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ABSTRACT 

We report detections of new exoplanets from a radial velocity (RV) survey of metal-rich FGK stars 
by using three telescopes. By optimizing our RV analysis method to long time-baseline observations, 
we have succeeded in detecting five new Jovian-planets around three metal-rich stars HD 1605, HD 
1666, and HD 67087 with the masses of 1.3Mg, 1.5M 0 , and 1.4M Q , respectively. A K1 subgiant 
star HD 1605 hosts two planetary companions with the minimum masses of M p sini = 0.96 Mjjjp 
and 3.5 Mjjjp in circular orbits with the planets’ periods P = 577.9 days and 2111 days, respectively. 

HD 1605 shows a significant linear trend in RVs. Such a system consisting of Jovian planets in 
circular orbits has rarely been found and thus HD 1605 should be an important example of a multi¬ 
planetary system that is likely unperturbed by planet-planet interactions. HD 1666 is a F7 main 
sequence star which hosts an eccentric and massive planet of M p smi = QAMjup in the orbit with 
a p = 0.94 AU and an eccentricity e = 0.63. Such an eccentric and massive planet can be explained as 
a result of planet-planet interactions among Jovian planets. While we have found the large residuals 
of rms = 35.6 m s~ 4 , the periodogram analysis does not support any additional periodicities. Finally, 

HD 67087 hosts two planets of M p sini = 3.1 Mjjjp and 4.9 Mjjjp in orbits with P = 352.2 days and 
2374 days, and e = 0.17 and 0.76, respectively. Although the current RVs do not lead to accurate 
determinations of its orbit and mass, HD 67087 c can be one of the most eccentric planets ever 
discovered in multiple systems. 

Subject headings: planetary systems — stars: individual(HD 1605, HD 1666, HD 67087) — techniques: 
radial velocities — techniques: spectroscopic 


1. INTRODUCTION 

Since 1995, precise radial-velocity (RV) searches for ex¬ 
oplanets have unveiled hundreds of planets around solar- 
type stars. The improvements of spectrographs and RV 
analysis techniques have enabled us to detect low-mass 
planets such as Earth-like planets. The current detection 
limit on planetary masses for RV searches is reaching 
the level of seve ral Earth-masses w i thin 0.1AU aroun d 
solar-type stars ((Howard et al.ll20Tol : iMavor et al.ll2011h . 
To date, masses of known exoplanets range from about 
one Earth-mass to tens of Jupiter-masses. Thanks to in¬ 
creasing the number of planet samples, statistical stud¬ 
ies become available. In particular, Jovian-planet sam¬ 
ples from 20-year-long surveys have been contributing to 
reveal their statistical properties over a wide range of 
orbital distance. The pictures of the detected exoplan- 
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ets have been contributing to constrain planet forma¬ 
tion theories taking account of evolutions of protoplan¬ 
etary disks and planet’s orbits such as planetary migra¬ 
tion an d planet-planet scattering (e.g. population syn¬ 
thesi s: llda fc Linll2004 llda et al.ll2013l : iMordasini et all 
2009), ho wever, statistical stud i es for the planet distribu¬ 
tion (e.g.lCumming et al.lll999ll20d8t lBowler et al.ll2010t 
iMortier et al.ll2012HAdibekvan et al .11 20131) are still open 
to debate, particularly on the correlation with metallici- 
ties of their host stars. Recently, near-infrared observa¬ 
tions of disk fractions in young clusters suggested short 
lifetim es o f p roto pla neta ry disks around low-metallicity 
stars (lYasui et all l2010lb Theoretical studies on disk 
evolution under X-ray photoevaporation supported that 
disk lifetimes may be positively correlated w ith the stel¬ 
lar metallicity (e.g. lErcolano &: Clarkelf2010f) . The disk 
metallicity also determines the total amount of build¬ 
ing materials for planets. These imply that a long-lived 
disk in a metal-rich environment may enhance the oc¬ 
currence rate of Jovian planets such as hot-Jupiters and 
Jupiter-analogs. While pr evious studies of the planet- 
metallicity correlation (e.g. IGonz alez 19 9*7t ISantos et all 

120031 iFischer fe Valentill2005H Johnson et al .1120 1 Pali have 

mainly focused on the occurrence rate of Jovian planets 
integrated over planetary mass and orbital period, it is 
also important to reveal dependence of the mass-period 
distribution of planets on stellar metallicity for an in- 
depth understanding of orbital migration and formation 
of giant planets in various environments. 

Since 2009, we have conducted RV surveys for metal- 
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rich solar-type stars mainly using the 188 cm telescope 
at Okayama Astrophysical Observatory (OAO) and the 
Subaru telescope. This planet-search program origi¬ 
nates from the N2K project, which was an interna¬ 
tional collaboration between American, Japanese and 
Chilean astronomers aiming at finding hot-Jupiters, es¬ 
pecially transiting ones, around metal-rich FGK stars us¬ 
ing large aperture telescope such as the Keck, Subaru , 
and Magellan (e.g. IFischer et al.ll2005HSato et. all [20051 

iLonez-Morales et all 120081 IGiguere et al.l l2012h . As for 
the Subaru telescope part, in addition to the detections 
of transiting HJs, continuous observations revealed Jo¬ 
vian planets from the including those in years-long or¬ 
bits (IFischer et al.ll2007t iPeek et aI1l2009t IGiguere et al.1 
l2012f ). The metal-rich nature of the N2K targets totally 
fits our purpose which aims at unveiling orbital distri¬ 
bution of giant planets around metal-rich stars over a 
wide range of orbital periods. In order to achieve this 
objective, long time baseline RV observations are of par¬ 
ticular importance in revealing a population of planets 
at distant orbits where Jupiter-analogs are thought to 
have formed via a core-accretion process. Continuous 
RV observations also lead to secure detec tions of mul¬ 
tiple systems and eccentric p lanets (e.g. iWright et all 
120091 : iWittenmver et all 12014) which will provide use¬ 
ful information about orbital evolutions such as an or¬ 
bital migration which occurs due to an interaction be¬ 
tween the disk-gas and planetesimals, and gravitational 
interactions among planetary /stellar compa nions (e.g. 
INagasawa et alJfefiost INagasawa fc Idall20lH) . 

In this paper, we report the detections of five new plan¬ 
ets around three metal-rich stars HD 1605, HD 1666, 
and HD 67087 from observations at telescopes of OAO, 
Subaru, and Keck. In Section [2j we present stellar char¬ 
acteristics and orbital solutions for the three stars. We 
also briefly describe the method for RV measurements 
including our improvements in RV stability and preci¬ 
sion for the Subaru telescope. In Section [3l we show our 
results of planet detections. We discuss and summarize 
the detections from a view point of planet formation and 
evolution in Section HI 

2. OBSERVATIONS 

2.1. Targets 

The stellar targets of our long-term planet-search pro¬ 
gram at OAO and the Subaru telescope are based on 
those of the N2K project. More than 14,000 main- 
sequence and s ubgiant stars were select ed from the Hip- 
parcos catalog dPerrvman fe ESA||1997|) as potential tar¬ 
gets of the N2K project, and their properties had been 
investigated such as B — V color magnitudes, JHK- 
band photometry from the 2MASS catalogue, stellar lu¬ 
minosity, parallaxes, proper motions, photometric varia¬ 
tions, and the existence of stellar companions. All the 
targets had been selected with following criteria: (1) 
0.4 < B — V < 1.2, (2) V < 10.5 and (3) stars closer than 
110 pc. The approximate 2,000 stars had been selected 
from the potential targets as the original N2K targets by 
setting metallicity criterion of [Fe/H] > +0.1 estimated 
from photometry. 

As the component of the N2K project, the Subaru tele¬ 
scope had observed a total of 635 stars out of the original 
N2K targets at least once as of 2009, and we adopted the 


635 stars as our primary targets for the subsequent ob¬ 
servations. Among them, we have intensively observed 
50 bright stars at OAO in order to confirm planet can¬ 
didates. Also, we have observed the 635 primary targets 
evenly at the Subaru telescope in order to mainly extend 
the time baseline of our samples for statistical studies. 

2.2. OAO/HIDES Observations 

For observations at OAO, we use the 188 cm tele¬ 
scope and the High Dispersion Echellc Spectrograph 
(HIDES; llzumiural 1 19991) installed at its coude focus. 
Until December 2010, we used the conventional slit aper¬ 
ture of HIDES (HIDES-S). The slit width was set to 
250^m(0".95), giving a spectral resolution R = 50000 by 
4 pixel sampling. A mosaic of three CCDs is capable of a 
simultaneous wavelength coverage from 3700 to 7500 A. 
For precise RV measur ements, we use d an iodine absorp¬ 
tion cell (U cell; IKambe et al.1 120021) . For V ~ 8 stars 
with an exposure time of 30 minutes, a typical signal- 
to-noise ratio SNR > 130 pix -1 at 5500 A, with which 
we can achieve the RV precision of 7 m s -1 . After Jan¬ 
uary 2011, we have been using the HIDES fiber-feeding 
system (HIDES-F) that has been newly developed for 
attaining better throughput (IKambe et all 12013 ). The 
HIDES-F uses an image-slicer as an entrance aperture 
of the spectrograph and its spectral resolution is fixed 
to R = 47, 000. We can typically obtain spectra with 
two times better efficiency (SNR > 180 pix -1 at 5500 A 
for V ~ 8 stars with a 30 min-exposure) compared to 
the previous slit-observations, with which the RV pre¬ 
cision can reach down to 3 ms -1 . We have been ob¬ 
serving several RV standard stars using both a slit and 
image-slicer, and monitoring these RV variations. As 
a Doppler-shift standard, we applied a newly obtained 
template-spectrum of HDS whose details are described in 
section m The calculated RVs did not show significant 
systematics between the slit-data and image-slicer-data 
within the typical rms (root-mean-square) of 3 ms -1 . 
Thus, we are able to consider that there is no significant 
RV systematics between the two settings of HIDES. 

Since the N2K project supposed planet searches by 10 
m-class telescopes, the apparent magnitude of the tar¬ 
gets tends to be fainter than 9 in V-magnitude. Stars 
with magnitudes of V > 9.5 are too faint for precise 
RV measurements by 2 m-class telescopes. Thus, we 
have selected 50 brightest stars that have magnitudes 
of 7.5 < V < 9.5 from our primary targets. Almost 
all of the 50 targets show significant RV variations of 
ctrv > 20 nr s^ 1 , but their periodicities have not been 
confirmed yet. Such RV variations are comparable to a 
RV amplitudes caused by a giant planet orbiting around 
a solar-mass star with a semi-major axis of 1 AU. We ob¬ 
tained stellar parameters of the 50 targets f rom the S pec¬ 
troscopy Made Easy (S ME; IValenti fe Piskunovl 119961 : 
IValenti fc Fischenl2005l) analysis and the Hipparcos cat¬ 
alogue. 

2.3. Subaru/HDS Observations 

From 2004 to 2010, the original N2K survey collected 
several RV data for each of our primary targets. Since 
2011, we have conducted follow-up observations for a 
part of them using the Subaru telescope. While observa¬ 
tions at OAO mainly focused on planet-host candidates 
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which showed significant RV variations. In order to ob¬ 
tain reliable statistical results, we need to gather RV data 
without selection biases. Thus, we collected RV data uni¬ 
formly from the primary targets as well as planet-host 
candidates, using the Subaru telescope. 

We ut ilize the High Dispersion Spectrograph (HDS; 

I Noguchi et al.ll200it) on the 8.2-m Subaru telescope. We 
mainly adopted the set-up of StdI2b, which simultane¬ 
ously covers a wavelength region from 3500 to 6100A by 
a mosaic of two CCDs. The slit width was set to 0".8 for 
the first 2 year-observations (2004-2005) and 0".6 for the 
later observations until 2011, giving a spectral resolution 
of R = 45000 and 60000, respectively. In 2012, we used a 
newly-developed image-slicer (IS#1) with the resolution 
of R = 110000 in order to improve the throughput of each 
exposure fsee lTaiitsu et al.l[2012t for details of the image- 
slicer). Our set-up allows us to obtain SNR ~ 150 pix -1 
at 5500 A with 30-300 sec exposure times for our typi¬ 
cal targets with V ~ 8.5. We used an I2 cell to provide 
a fiduci al wavelength reference for precise R V measure¬ 
ments (iKambe et all 120021 iSato et all I2002H . Although 
the I2 cell was originally installed just behind the en¬ 
trance slit of the spectrograph, it was moved in front of 
it in 2006. 

We basi cally used an analysis code described in 
ISato et ahl (20021 1 for precise RV measurements which 
generates the stellar template from star+12 spectra. Stel¬ 
lar templates obtained in the early phase of the N2K 
survey sometimes cause poor RV accuracy for years-long 
time-baseline observations with HDS. We also conducted 
additional observations of the N2K targets in order to re¬ 
fine stellar templates. The observational setting required 
for the stellar templates is almost the same as that for 
RV measurements, except for using extremely high reso¬ 
lution of R = 160000 (the slit width of 0".2 without an 
I2 cell). The point-spread-functions (i.e. instrumental- 
profiles) are estimated from an ^-superposed flat spec¬ 
trum that was observed at the nearest time before/after 
the exposure. Using the newly-obtained stellar tem¬ 
plates, we can achieve the rms of RVs lower than 3ms -1 
with about seven-year-long stability. We optimized our 
analysis pipeline for this refined template, and applied 
it to RV measurements for HIDES data as a common 
reference of Doppler-shifts. 

In order to check the long-term RV stability using the 
whole data from HIDES and HDS, we have been observ¬ 
ing several RV sta ndard stars and know n planet hosts 
such as HD 38801 (lHarakawa et all [20101 1 using the two 
spectrographs. According to our analysis for ”slit-IS” 
combined set of data from HDS, it has been shown that 
there is no significant systematic between slit data and 
image-slicer data. We have also confirmed that the 12-cell 
position against the entrance of HDS does not influence 
RV systematics significantly. After these checks, we have 
performed the RV measurements for all the HIDES and 
HDS data and confirmed that the rms of RVs for the 
standard stars have been achieved 4-5 m s _1 . Thus, we 
can consider that no significant RV-offset between the 
two sites was detected with a mixture of slit-data and 
image-slicer-data. 

2.4. Keck/HIRES Observations 


Between 2006 January and 2010 September, we also 
used t he Keck telescope with HIRES spectrograph 
(Vogt et al.l 11994T) to follow up some promising candi¬ 
dates of planet-host stars. The Keck observations were 
obtained using an I2 cell and typically with a slit width 
of 0”.86 corresponding to a spectral resolution of R = 
65000. The exposure time is optimized to achieve the 
typical signal to noise ratio of 150-240 per pixel. The 
mean RV precision of 1.5m s -1 was yielded based on th e 
RV analysis method pro posed bv lMarcv fc Butleul (jl992), 
and lButler et al.l ([1996). Since the RVs of HIRES were 
calculated individually from those of HIDES and HDS, 
the RV offsets between the two datasets are fitted as free 
parameters. 

3. STELLAR CHARACTERISTICS AND ORBITAL 
SOLUTIONS 

3.1. HD 1605 

HD 1605 (HIP1640) is a K1 subgiant star with V = 
7.52 and a color index of B — V = 0.96. The Hip- 
parcos parallax of 11.82 mas corresponds to a distance 
of 84.6 pc, i.e., the absolute V magnitude of M v = 
2.88. The SME analysis yields the stellar parameters 
of T eff = 4757 ± 50 K, vsmi = 0.54 ± 0.5 km s' 1 
[Fe/H]= +0.21 ±0.05, and log g = 3.40±0.08. Adopting 
the determined T e g , B — V , My, and bolometric collec¬ 
tion of B.C. = -0.432 (|Flowerill99fifl . we have estimated 
the bolometric luminosity to be L = 6.6 ± O.76L 0 , and 
stellar radius of i?* = 3.8 ± O.4i? 0 based on the Stefan- 
Boltzmann law. We estimated the stellar mass from the 
be st-fit i nterpolation of stellar evolution tracks calculated 
bv lGirardi et, al.l ([2002H to be M* = 1.31 ± O.IIMq with 
its age of 4.59 ± 1.37 Gyr. In fact, [Fe/H]> +0.2 is out 
of range of availabl e data of stellar evolution models by 
iGirardi et ahl (120021) . Thus, we have extrapolated the ex¬ 
ist ing_tracks_for models with a higher stellar metallicity 
(see lOmiva et al.ll2012l for the treatment of our extrap¬ 
olation method). All the stellar parameters are listed in 
Table [Q 

We made the first observations for HD 1605 at the 
Subaru telescope and detected large RV variations. After 
that, the star has been monitored using HDS, HIDES 
and HIRES for about 9 years (from December 2004 to 
August 2013). We excluded data with low SNR under 
bad weather conditions, and then obtained 75 RV data, 
14 and 61 of which were observed by HDS and HIDES, 
respectively. Combining 17 RV data from HIRES, we 
have obtained 92 RV data in total for HD 1605. The 
observation dates, RVs, internal errors, and observation 
sites are listed in table [2] 

We determined parameter values of the Keplerian or¬ 
bit by minimizing th e x 2 statistic. We adopted the 
AMOEBA algorithm (|Nelder fc Meadl 1 1 9651 [Press et ahl 
1986) which is well-suited to the high-dimensional pa¬ 
rameter spaces. The definition of x' 2 is written as 

Nsite A+obs , \ 2 

2 -L I ^ ij ,obs calc T7i \ 

-^total -^dim \ error / 

i j \ ji / 

(i) 

where N sitei 

,obS5 ^ ij , obs? & ij , error •> ^ ij , calc5 

and 7 i are the number of observational sites, the num¬ 
ber of observations on the *-th site, the number of to¬ 
tal observations, the degree of freedom on the fitting, 
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the j- th observed RV and observational uncertainty from 
each site, the calculated RVs from the model-orbit cor¬ 
responding to j- th observation time on i-th site, and RV 
offsets on the *-th site, respectively. In fact, the number 
of observational sites can be written as the number of 
stellar templates. A quantity of the RV offset depends 
on the time when the stellar template was observed. In 
this work, we used not more than two stellar templates 
in total that had been obtained from HDS and HIRES. 
Therefore, the appropriate number of site is V g ^ e = 1 
or 2. Each template was applied to RV measurements 
of the HDS-HIDES combined data set, and HIRES data 
set, respectively. 

In order to identify the stellar chromospheric activity, 
we graphically checked the core emission strength of Ca 
II HK lines for each star. As shown in Figure [TJ all 
of the stars are basically inactive. We assumed an in¬ 
trinsic RV noise due to the chromospheric activity of a 
central star (e.g. starspots and plage) and its surface 
pulsation, the so-called “stellar jitter”. to b e 6 ms -1 
for ordinary inactive-subgiants (IWrightl [200511 . which is 
added in quadrature (<7g rror = fj^ bs + cr? itter ) to the for¬ 
mal velocity uncertainties for HD 1605 when we fit a 
Keplerian orbit to the R V data. At first, we perfo rmed 
the periodogram analysis dbombll 197(1 lScargldfi982H and 
found significant periodicity at ~ 2000 days. After de¬ 
trending the 2000 days signal, an another periodic signal 
appeared at several hundred days period. Although we 
considered a two-planet model for the Keplerian fitting, 
its residuals showed a linear slope over the entire obser¬ 
vations. Therefore, we derived the best-fit solution of a 
two-planet model with a linear trend. The rms to the 
fit is reduced from 10 m s -1 to 6.4 m s _1 . The reduced 
chi-square has been improved from 1.9 (2 planet-model) 
to 0.94 (2 planet-model + linear trend), indicating that 
the latter model should be more appropriate within our 
achievable accuracy. We calculated the errors in orbital 
parameters at a 68.27% confidence interval (from bottom 
15.87% to top 84.13%) of the distributions via the boot¬ 
strap Monte-Carlo approach. The RV residuals to the 
best-fit Keplerian curve are scrambled and added back 
to the original measurements, and then, we re-fit and 
derive new parameters. We achieved la confidence lev¬ 
els of each parameter by 10000 times iterations. We also 
applied the same technique to two other systems reported 
in the subsequent sections. 

Figure [2] show the RV variations and the RV curve 
based on best-fit values. The best-fit orbital solutions 
and their errors are listed in tablc[5] The orbital solutions 
of planet “b” are P = 577 . 9 + 4.9 days and the RV semi¬ 
amplitude A'i = 19.8l 0 'Jm s_1 ’ corresponding to the 
minimum planetary mass M p sin i = 0.96 Mjjjp at the 
semimajor axis a p = 1.48 AU, and those of outer planet 
“c” are P = 2111 +32 days and K\ = 46.5j)i‘Jm s_1 ) 
respectively, corresponding to M p sin i = 3A8Mjjjp at 
a p = 3.52 AU. The remaining RV trend shows the linear 
RV variation of —4.92 m s _1 yr -1 . The orbital period ra¬ 
tio of the planets is not any integer even considering their 
ler uncertainties. This suggests that two Jovian planets 
are not trapped into any mean motion resonance. Both 
of the two planets have almost circular orbits (e < 0 . 1 ), 
and thus we cannot strictly constrain their eccentricities 
and arguments of periapsis. 


3.2. HD 1666 

HD 1666 (HIP1666) is listed as a F7V star with a visual 
magnitude of V = 8.17, and a color index B — V = 0.53 
in the Hipparcos catalogue. The Hipparcos parallax of 
9.22 mas corresponds to a distance of 108.5 pc, and the 
absolute visual magnitude is estimated to My = 2.76. 
We performed the SME analysis in order to determine 
th e stellar paramete r s. Fo llowing the method described 
in IValenti fc Fischeil (12005D . we derived an effective tem¬ 
perature T e ff= 6317 ±44K, rotational velocity wsinz = 
5.6 ± 0.5km s _1 , log g = 4.06 ± 0.06, and metallicity 
[Fe/H]= +0.37 ± 0.04 for the star. The bolometric lu¬ 
minosity of L* = 5.37 ± 1.18 -Lq was estimated from the 
My and a bolometric correction of B.C. = —0.016. We 
derived a stellar radius of i?* = 1.93 ± O.49i ? 0 from 
the Stefan-Boltzmann law. We performed the isochrone 
analysis as mentioned in Section 13.11 based on the stellar 
metallicity, bolometric luminosity and effective tempera¬ 
ture. Using the spectroscopic parameters, we found the 
mass of 1.50 ± O.O7M 0 and the age of 1.76 ± 0.20 Gyr for 
HD 1666. 

After the detection of large RV variations for HD 1666 
by HDS, we have monitored the star for about 9 years 
(from December 2004 to August 2013) using the three- 
spectrographs. We excluded data with extremely-low 
SNRs due to bad weather conditions and obtained 99 
RV measurements in total, 11, 21, 46, and 21 of which 
were observed by HDS, HIDES-S, HIDES-F, and HIRES, 
respectively. Table [3] summarizes the observation dates, 
RVs, internal errors, and observatory sites. 

The observed RVs showed complex patterns, as seen 
in Figure [3] We applied a Lomb-Scargle periodogram 
to the measured RVs of HD 1666 for the determination 
of the initial guess of orbital parameters and the result 
is shown in Figure |U We found a ~ 270 days period¬ 
icity with a confidence level higher than 99.9% by cal¬ 
culating the False-Alarm-Probability (FAP) defined by 
iBaluevI (|2008f ). In order to derive the Keplerian orbital 
motion for HD 1666, we adopted the AMOEBA method 
as described in the previous section. This star should be 
chromospherically inactive bec ause of the small activity- 
index value of Shk = 0.15 (llsaacson fc Fischeil 1201 04 
and the low-strength of the cores of Ca II H K lines 
(Figure [T]). According to an empirical eq uatio n o f ste llar 
jitters, as a function of Shk given bv llsaacson fc Fischerl 
(IMP , we estimated the RV jitter of HD 1666 to be 
2.4 m s -1 . However, as seen in Figure^ the RV-residuals 
of the orbital fitting remain still large variation in rms 
of > 30 m s -1 , compared to the typical observational 
uncertainty of 12 m s -1 . In general, stellar RV varia¬ 
tions are caused by both the Keplerian motion and stel¬ 
lar jitter. The 104 data from the Hipparcos photometry 
(|Perrvman fc E SA 1997) show a photometric variability 
of a ~ 13 mmag for an inactive HD 1666. However, 
this value is comparable to the typical observational er¬ 
ror of 15 mmag, indicating that there is no significant 
variation photometrically for HD 1666. This implies that 
there may exist additional companions that are respon¬ 
sible for the excess RV variations. While a period search 
for the residuals to the fit detected a slightly strong signal 
at about ten days (bottom panel of Figure [4]), the FAP 
value of 3.5% does not strongly support its significance. 
Indeed, although we tried an additional Keplerian-orbit 
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fitting to the residuals with an initial guess of ten days 
period, neither the reduced chi-square nor the rms was 
improved. More precise and additional data are required 
to interpret these large RV-residuals. Then, for the 
present, we regarded these large residuals as the stel¬ 
lar jitter, (Jjitter = 34.6 m s -1 , which makes the reduced 
chi-square unity. Finally, we obtained the orbital pe¬ 
riod of P = 270.0+o)g days, the velocity semi-amplitude 
K\ = 199.4^5’ 8 m s _1 , and the eccentricity e = 0.63+0+2 
as the best-fit orbital parameters for HD 1666 b. With 
the stellar mass of 1.50 ± 0.07Mg, we found that HD 
1666 b has M p sini = 6A3Mjpp at the semi-major axis 
of a p = 0.94 AU. The best-fit parameters and their errors 
are listed in Table [5] 

3.3. HD 67087 

HD 67087 (HIP39767) is a F-type main-sequence star 
with V = 8.05 and a color index of B — V = 0.53. 
The Hipaarcos parallax of 11.26 mas corresponds to a 
distance of 88.8 pc, i.e., the absolute V magnitude of 
My = 3.31. The SME analysis yields the stellar param¬ 
eters of T e g = 6330 ± 47 K, v sin i = 9.7 ± 0.5 km s -1 , 
[Fe/H]= +0.25 ± 0.04, and log g = 4.19 ± 0.06. Using 
T e g, B — V, My, and B.C. = —0.010, we have esti¬ 
mated the bolometric luminosity and stellar radius to be 
L = 3.47 ± 0.75 Lq and A* = 1.55 ± O.32i? 0 from the 
Stefan-Boltzmann law. The stellar mass and age were 
estimated to be M* = 1.36 ± O.O4M 0 and 1.45 ± 0.51 
Gyr from stellar evolution tracks, respectively, (see also 
Section 13.11) . 

We have monitored HD 67087 using HDS and HIDES 
for about 9 years (from December 2004 to May 2014). We 
obtained 51 RV data in total, 12 and 39 of which were 
observed by HDS and HIDES-F excluding poor observa¬ 
tions, respectively. The observation dates, RVs, internal 
errors, and observatory sites are listed in Table [4j The 
observed RVs show a clear significant periodicity in ad¬ 
dition to that of about one-year in Figure[5l Assuming 
the stellar jitter of 4 m s _1 based on the low emission 
strength in the core of Ca II HK lines (Figured]), we fit¬ 
ted two-Keplerian orbits to the measured RVs. However, 
we cannot obtain the RV semi-amplitude and the min¬ 
imum mass of a planet, via the bootstrap method due 
to lack of RV peaks for the outer orbit. Hence, we ap - 
plied the MCMC procedure following iSato et all ( 2013 ) 
to the RVs and determined two-Keplerian orbits from the 
posterior distribution functions. 

We considered the medians at the la confidence level 
as the best-fit values. The rms to the fitting was 
11.8 ms 1 , and the square root of reduced \ 2 of \fx^, = 
0.91, indicating that our orbital solutions are prefer¬ 
able within observational errors. With the stellar mass 
of M* = 1.36 ± O.O4M 0 , P = 352.27±};£, and K\ = 
73 . 6 l^m s-\ we found a p = 1.08 AU and M p sini = 
3.06Mjup for the inner planet HD 67087 b. The outer 
planet HD 67087 c has the orbital period P = 2374l^g 
days in an eccentric orbit of e = 0 . 76 ^ 0)24 and K\ = 
OS.Sl^g 4 m s-\ and hence, M p sini = 4.85tg ° 6 iMjup 

at a p = 3.86+0 37 AU. Figure [5] shows observations and 
best-fit orbital solutions. The orbital solutions and their 
errors are summarized in Table [5] 

Although the phase-folded RV variation of HD 67087 c 


appears a simple linear RV trend, the actual time varia¬ 
tion in RVs repeats this trend more than one within the 
duration of observations. We tried the Keplerian fitting 
with one-planet + linear trend model and obtained a re¬ 
duced chi-square of xt = 3-4. Therefore, the two-planet 
model may be preferred for this system, considering the 
reduced chi-square of xl = 0.81. 

4. DISCUSSION 

In this paper, we have reported five newly-discovered 
Jovian planets around metal-rich stars from our long¬ 
term planet-search at the Subaru telescope, the OAO 
188 cm telescope and the Keck telescope. Two multi- 
Jovian-planet systems of HD 1605 and HD 67087, and 
the eccentric massive planetary system of HD 1666 have 
been identified. All of the host stars have high metal- 
licity ([Fe/H]> +0.2) and relatively large masses of 
M* > 1.3M 0 . 

HD 1605 is a subgiant star and its mass is estimated 
to be 1.3 Mq. We have confirmed the two Jovian-mass 
planets on the hierarchical system with a linear RV trend. 
According to the SIMBAD database, no stellar compan¬ 
ions has been reported within 1.1 arcsec (<100 AU) from 
HD 1605. The linear trend of RV that we found suggests 
the presence of an unknown and maybe low-mass com¬ 
panion. Given that the acceleration of RVs is constant, 
we can give an order-of-magnitude relation between the 
linear trend and the properties of the outer companion 
(IWinn et all 1 20091) . 

„ 1 = (0.028AU" 2 , (2) 

a c Lr 

where M c is the companion mass, i c is the orbital inclina¬ 
tion relative to the line of sight, a c is the orbital distance, 
and 7 is the acceleration of RVs. The HD 1605 system 
is also a rare planetary system in terms of their low- 
eccentricity orbits. There are 44 systems which hold mul¬ 
tiple Jovian-planets with masses of M n smi > 0 .3 Mjpp 
(the Exoplanets Orbit Database: iWright et al.l 1201 ill . 
and only four of them are known to have multiple Jo¬ 
vian planets with low-eccentricities (e < 0 . 1 ) in the 
absence of h ot Jupiters (47 UMa, 55 Cnc, p Ara and 
HD 159243: iGregorv &; Fischer] 120101 : lEndl et all 120121 
iPeoe et alll2007HMoutou et al.ll2014 respectively). 

The two-planet system of HD 67087 is also a hierar¬ 
chical system. HD 67087 has the mass of 1.36M 0 and 
metallicity of [Fe/H] = +0.25. We have found two super- 
Jupiter-mass planets at 1.08 AU and 3.86 AU around HD 
67087. Interestingly, HD 67087 c may be one of the most 
eccentric Jovian-planets among those with a p > 3 AU 
and also among multiple Jovian-planet-systems, in con¬ 
trast to HD 1605 c (Figure ©• Although Jovian planets 
similar to HD 67087 c a re seen in o ther systems such 
as HD 181433 d (IBouchv et al 1^ 0091 and HD 217107 c 
(|Vogt et al.l 120051 : 1 Wright et alJl2009h . hot-Jupiters were 
found in their inner orbits. Thus, the origin of those ec¬ 
centric outer planets can be explained by planet-planet 
interactions. On the other hand, in the case of HD 67087 
which has no hot-Jupiter, the existence of an eccentric 
outer planet around HD 67087 may be interpreted by 
the excitation of its eccentricity due to secular pertur¬ 
bations from an unseen substellar co mpanion in distant 
orbit, the so-called Kozai mechanism (|Kozailfl962f ). Al¬ 
though the lack of observations does not lead to the ac- 





















6 


Harakawa et al. 


curate determination of the orbit of HD 67087 c, the two 
planets in this system have mutually crossing orbits at 
present. Hence, the dynamical stability of HD 67087 
system should be investigated with additional RVs in 
the future work. This allows us to constrain the mass 
and the eccentricity of HD 67087 c more strictly. Al¬ 
though we did not find any long-term RV trend due to 
an additional companion from the residuals to our best- 
fit model, a relatively-young HD 67087 still becomes a 
promising target for direct-imaging surveys. If the exis¬ 
tence of an outer companion is confirmed, HD 67087 c 
should be a strong evidence of orbital evolution via Kozai 
mechanism. 

The F7 dwarf HD 1666 has a relatively-large mass of 
1.5 -Mq and a very high metallicity of [Fe/H] = +0.37. 
HD 1666 hosts a massive, eccentric planet HD 1666 b. 
One of the most remarkable features of HD 1666 sys¬ 
tem is the mass of a host star. It is generally known 
that early-type main-sequence stars which have masses 
higher than that of the Sun such as A-dwarfs (typi¬ 
cally 1.5 Mq— 2Mg) are not suitable for RV measure¬ 
ments because of Doppler-broadening of absorption lines 
due to their rapid rotation and the deficiency of absorp¬ 
tion lines due to high effective temperatures. RV sur¬ 
veys rarely detected sh ort-period plan ets aro und evolved 
high-mass stars fe.g. I Johnson et ahl 120103) . However, 
HD 1666 has a cooler effective temperature and lower 
rotational velocity. A high metallicity in the stellar at¬ 
mosphere causes a larger stellar radius at a given lumi¬ 
nosity and stellar mass. A metal-rich star like HD 1666 
with [Fe/H] = +0.37 tends to have a lower effe ctive 
temperature (e.g. lGirardi et al.ll2000HYi et al.ll200llf . As 
seen in Figure 17 of Valenti & Fisclien~ it2005ll . T e g is de¬ 
creased by about 500K, when [Fe/H] increases by 0.5 
dex. Our target selection for a metal-rich star based on 
the color-index (typically B—V) leads to selecting a high- 
mass star, and vice versa. Figure [7] shows the planetary 
mass distribution as a function of stellar mass. All of 
the Jovian planets in this figure show those around main 
sequence stars (logg > 4.0) detected by the RV method. 
Figure [7] indicates that HD 1666 is one of the most mas¬ 
sive main-sequence stars harboring planetary compan¬ 
ions. Therefore, metal-rich and massive main-sequence 
stars such as HD 1666 can be suitable targets not only to 
unveil a correlation between planetary properties and the 
stellar metallicity, but also to provide information about 
the dependence of planet population on stellar mass. 

Another important feature of HD 1666 system is the 
eccentric planet of HD 1666 b. As mentioned before, such 
an eccentric planet can be interpreted as a by-product of 
gravitational interactions among planets and/or stellar 
companions, although we have not yet found any inner 
or o uter co mpanion. The detection-limit analysis follow¬ 
ing [Lagrange/erid] (120091) is shown at Figure [S] We here, 
briefly describe about this analysis method: (1) for each 
(mass; period) grid, generate 1000 synthetic RV-data- 
sets along the actual observation times assuming circu¬ 
lar orbit and different times of passage at periastron; (2) 
add a random noise corresponding real observational er¬ 
rors; (3) compare the averaged rms of virtual data-sets 
and that of real residual data. If the rms of the real 
residuals is lower than the average rms of virtual data 
sets, the given planet can be detected. The detection 
probability can be calculated as a level of rms excess 


compared to that of the real data. Figure [8] shows that 
while the stellar jitter (emitter = 35 nr s” 1 ) may prevent 
us from detecting low-mass planets, it cannot be an ob¬ 
stacle for the detection of a hot-Jupiter with the mass 
of MpSini > 3 Mjjjp with P < 10 days at a 3er confi¬ 
dence level. As seen in Figure [2 the non-detection for 
short periodicity from the residuals indicates that (i) if 
HD 1666 b experienced planet-planet scattering, an scat¬ 
tered planet is likely to be outside HD 1666 b or (ii) the 
eccentricity of HD 1666 b was pumped up by an outer 
companion via the Kozai mechanism. Indeed, the RV 
amplitude induced by a distant substellar companion can 
be smaller than this jitter. The future direct-imaging sur¬ 
vey or correction for stellar jitters with careful analyses of 
activ i ty-induced RVs (e.g.lPont■ et al.ll2011k iBoisse et al.1 
I2012t IHatzesll2013l : ISantos et al.ll2014ll helps us verify the 
existence of the additional companion(s). 
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Wavelength (A) 


Figure 1 . Spectra centered on the Ca II H line for the three planet-host stars. All of the stars show no significant emission in the core of 
lines indicating chromospherically quiet. An appropriate offset is added to each spectrum for clarity. 
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4000 5000 6000 


JD (-2450000) 



Orbital Phase 

Figure 2. Top panel: Observed RV variations and the best-fit orbital solution of HD 1605 with its residuals to the fit. Light-blue points, 
green points, blue points, and pink points represent data from the Subaru/HDS, OAO/HIDES-S, OAO/HIDES-F, and the Keck/HIRES, 
respectively. The best-fit stellar RV motion has been shown as red solid line. The error bars are shown as the summations of quadrature 
of observational errors and the stellar jitter of ajitter = 6 m s —1 . The dashed line shows a linear trend that indicates the existence of 
uncovered (sub)stellar or planetary companion(s). The offsets of 6.7 m s —1 and 0.1 m s —1 were introduced to HDS and HIDES data, and 
HIRES data, respectively. Bottom panel: Phase-folded RV variations for HD 1605. 
















10 


Harakawa et al. 



4000 5000 6000 

JD (-2450000) 



Orbital Phase 

Figure 3. Top panel: Observed RV variations and the best-fit orbital solution of HD 1666 with its residuals to the Best-fit-model. Error 
bars are drawn as the summations of quadrature of observation uncertainties and assumed stellar jitter of Oji tter = 34.6 ms -1 . The RV 
offsets of —7.4 m s —1 and 54.3 ms -1 were added to HDS and HIDES data, and HIRES data, respectively. The color legends are the same 
as those of Figure [2] Bottom panel: Phase-folded RV variations of HD 1666. 
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Figure 4. Upper panel: the Lomb-Scargl e periodogram for RV measurements on HD 1666. FAP on the maximum peak 270 days) 
was calculated based on the definition from [Baluev! ^20081 ). Lower panel: same for the residuals to the Keplerian fit. 
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4000 5000 6000 

JD (-2450000) 



Orbital Phase 

Figure 5. Top panel: Observed RV variations and the best-fit orbital solution of HD 67087 with its residuals to the best-fit model. The 
color legends are the same as Figure [2] Bottom panel: Phased RV variations of HD 67087. 
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semimajor-axis [AU] 

Figure 6. Planet distribution in a p -e diagram. The planets with M p smi > 0.1 Mjjjp in multiple systems are shown. Only the error bars 
of HD 67087 planets are drawn for clarity. —http://exoplanets.org 
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M* [M sun ] 

Figure 7. Distribution of minimum masses of exoplanets and masses of their host stars. Only Jovian planets around main-sequence stars 
(Mpsini > O.IMjup, log# > 4.0) are plotted —http://exoplanets.org Planets around HD 1666 and HD 67087 are shown with filled red 
circle and filled orange circle , respectively. 
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Figure 8. Detection probability map calculated from the RV residuals of HD 1666. The red solid line and the black dashed line represent 
3cr(99.73%) and lcr(68.27%) detection probability, respectively. 
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Table 1 

Respective Stellar Parameters 


Parameter 

HD 1605 

HD 1666 

HD 67087 

V 

7.52 

8.17 

8.05 

M v 

2.88 

2.76 

3.31 

B — V 

0.96 

0.53 

0.53 

Spectral type 

K1 IV 

F7 V 

F8 V 

Distance (pc) 

84.6 

110.6 

88.8 

Teff (K) 

4757 ± 50 

6317 ± 44 

6330 ± 47 

log a {cgi) 

3.40 ± 0.07 

4.06 ± 0.06 

4.19 + 0.06 

[Fe/H] 

+0.21 ± 0.05 

+0.37 + 0.04 

+0.25 ± 0.04 

v sini (km s 1 ) 

0.54 ± 0.5 

5.6 ±0.5 

9.7 ±0.5 

M, (M 0 ) 

1.31 + 0.11 

1.50 + 0.07 

1.36 + 0.04 

R*(Rq) 

3.8 ±0.4 

1.93 + 0.49 

1.55 + 0.32 

L*{Lq) 

6.6 + 0.76 

5.37 + 1.18 

3.47 + 0.75 

B.C. 

-0.432 

-0.016 

-0.010 

age (Gyr) 

4.59 + 1.37 

1.76 + 0.20 

1.45 ± 0.51 
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Table 2 

RVs for HD 1605. 


Observation Date 
(JD -2450000) 

RV 

(ms -1 ) 

Uncertainties 

(ms -1 ) 

Observatory 

3242.026396 

55.28 

3.85 

Subaru 

3360.824924 

74.05 

3.47 

Subaru 

3361.734252 

68.41 

3.59 

Subaru 

3401.755798 

60.49 

3.69 

Subaru 

3693.825842 

6.87 

3.37 

Subaru 

3726.799867 

8.19 

3.60 

Subaru 

3727.684276 

5.18 

3.59 

Subaru 

3727.742699 

4.19 

3.36 

Subaru 

3727.823795 

-5.14 

3.56 

Subaru 

3746.791000 

17.27 

0.99 

Keck 

3747.788000 

2.04 

1.09 

Keck 

3748.716000 

12.47 

1.02 

Keck 

3749.719000 

2.58 

1.02 

Keck 

3750.706000 

4.03 

1.04 

Keck 

3775.734000 

17.89 

1.03 

Keck 

3776.712000 

6.55 

1.06 

Keck 

4427.814000 

-48.43 

1.07 

Keck 

4428.839000 

-50.16 

1.09 

Keck 

4635.073000 

-28.65 

1.01 

Keck 

4638.113000 

-33.57 

0.99 

Keck 

4641.098000 

-33.98 

1.10 

Keck 

4743.946902 

-40.84 

3.06 

Subaru 

4786.941025 

-51.52 

3.48 

Subaru 

5044.113000 

12.41 

1.06 

Keck 

5056.151530 

32.23 

4.63 

OAO/S 

5082.147000 

19.70 

1.01 

Keck 

5099.086084 

32.55 

3.62 

OAO/S 

5130.997159 

46.44 

4.90 

OAO/S 

5159.174942 

40.19 

4.61 

OAO/S 

5171.854000 

36.05 

1.07 

Keck 

5184.968137 

44.59 

3.66 

OAO/S 

5186.975068 

37.95 

4.70 

OAO/S 

5187.791000 

40.23 

1.07 

Keck 

5234.905796 

24.64 

15.17 

OAO/S 

5379.092000 

3.62 

1.09 

Keck 

5398.230411 

15.43 

4.54 

OAO/S 

5401.227137 

14.15 

7.79 

OAO/S 

5402.233880 

7.37 

4.29 

OAO/S 

5404.234587 

7.69 

4.67 

OAO/S 

5409.228767 

15.62 

3.72 

OAO/S 

5410.204903 

14.99 

4.46 

OAO/S 

5428.071809 

17.37 

3.22 

OAO/F 

5429.125070 

17.21 

3.22 

OAO/F 

5430.151899 

20.53 

3.14 

OAO/F 

5431.302805 

13.18 

3.05 

OAO/F 

5432.148032 

17.45 

3.20 

OAO/F 

5433.162621 

15.67 

3.53 

OAO/F 

5434.218731 

13.32 

3.45 

OAO/F 

5435.173909 

13.84 

3.06 

OAO/F 

5441.181714 

3.31 

3.11 

OAO/S 

5445.266316 

4.96 

4.13 

OAO/S 

5468.117129 

0.50 

4.45 

OAO/S 

5474.112849 

3.61 

3.55 

OAO/S 

5476.138541 

1.69 

3.52 

OAO/S 

5480.204613 

-2.43 

4.52 

OAO/S 

5502.094646 

5.50 

4.11 

OAO/S 

5506.089555 

1.85 

3.29 

OAO/S 

5525.033422 

3.43 

3.85 

OAO/S 

5545.010623 

7.29 

5.98 

OAO/S 

5546.990257 

18.82 

5.16 

OAO/S 

5570.972884 

27.02 

4.78 

OAO/F 

5596.902826 

25.46 

2.98 

OAO/F 

5618.910983 

25.10 

3.27 

OAO/F 

5696.303752 

47.87 

4.20 

OAO/F 

5697.295470 

39.86 

6.77 

OAO/F 

5700.291167 

25.83 

3.38 

OAO/F 

5701.292900 

30.94 

3.53 

OAO/F 

5757.160078 

30.11 

3.85 

OAO/S 

5758.277612 

11.75 

4.92 

OAO/S 

5759.162776 

22.78 

4.76 

OAO/S 

5772.170736 

29.60 

7.17 

OAO/F 

5772.212740 

19.76 

2.98 

OAO/F 

5781.135733 

5.52 

3.63 

Subaru 

5812.030734 

8.19 

2.79 

OAO/F 
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Table 2 — Continued 


Observation Date 
(JD -2450000) 

RV 

(ms -1 ) 

Uncertainties 

(ms -1 ) 

Observatory 

5813.038857 

8.60 

2.64 

OAO/F 

5814.041536 

11.69 

2.97 

OAO/F 

5817.074067 

9.19 

2.99 

OAO/F 

5842.081940 

-4.83 

2.94 

OAO/F 

5843.009333 

-2.09 

3.07 

OAO/F 

5910.995010 

-32.82 

3.77 

OAO/F 

5931.934054 

-35.43 

2.65 

OAO/F 

6132.215925 

-50.30 

5.02 

OAO/F 

6135.173177 

-40.70 

3.73 

OAO/F 

6138.167379 

-43.62 

4.32 

OAO/F 

6148.110402 

-40.38 

3.13 

Subaru/IS#l 

6149.102357 

-48.86 

4.70 

Subaru/IS#l 

6163.122505 

-38.39 

5.08 

OAO/F 

6240.975244 

-41.49 

3.91 

OAO/F 

6279.939767 

-42.60 

3.49 

OAO/F 

6514.129844 

-98.55 

3.97 

OAO/F 

6517.149862 

-96.22 

3.50 

OAO/F 

6523.128885 

-100.38 

3.62 

OAO/F 


Note. — OAO/S and OAO/F represent data observed from OAO 188 cm telescope with HIDES-Slit mode and 
HIDES-Fiber feed mode, respectively. The RV data from Subaru/HDS using IS#1 are presented as Subaru/IS#l. 
All RV values are subtracted by appropriate 7 value listed in Table [5] Subsequent RV tables follow this manner. 
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Table 3 

RVs for HD 1666. 


Observation Date 
(JD -2450000) 

RV 

(ms -1 ) 

Uncertainties 

(ms -1 ) 

Observatory 

3360.701854 

34.89 

6.71 

Subaru 

3361.693030 

11.01 

7.51 

Subaru 

3362.692462 

26.71 

6.80 

Subaru 

3401.709260 

-80.49 

9.41 

Subaru 

4078.817044 

281.88 

6.86 

Subaru 

4080.798905 

257.23 

7.26 

Subaru 

4081.784028 

237.31 

5.83 

Subaru 

4743.990920 

-80.56 

8.37 

Subaru 

5099.134235 

-139.93 

11.31 

OAO/S 

5131.047738 

-125.16 

16.27 

OAO/S 

5184.894363 

159.10 

12.13 

OAO/S 

5186.903685 

168.04 

11.78 

OAO/S 

5396.121716 

-191.64 

1.72 

Keck 

5397.063226 

-201.03 

1.64 

Keck 

5400.127388 

-170.83 

1.76 

Keck 

5401.051892 

-150.05 

2.00 

Keck 

5401.274669 

-126.25 

30.18 

OAO/S 

5402.283521 

-141.41 

14.87 

OAO/S 

5403.104754 

-108.01 

1.58 

Keck 

5404.282251 

52.41 

20.70 

OAO/S 

5405.092260 

-19.88 

1.73 

Keck 

5406.088778 

-53.29 

1.67 

Keck 

5407.107137 

-68.75 

1.62 

Keck 

5410.287100 

-11.85 

21.87 

OAO/S 

5411.112017 

-8.15 

1.56 

Keck 

5412.006057 

-20.98 

1.62 

Keck 

5413.071693 

19.40 

1.58 

Keck 

5414.079847 

91.11 

1.64 

Keck 

5426.081973 

295.80 

1.95 

Keck 

5428.102646 

286.88 

1.85 

Keck 

5429.100236 

272.53 

1.81 

Keck 

5429.220400 

242.83 

14.74 

OAO/F 

5430.199245 

219.08 

14.58 

OAO/F 

5431.210050 

211.42 

12.46 

OAO/F 

5432.241269 

177.79 

14.56 

OAO/F 

5433.243331 

129.94 

13.56 

OAO/F 

5433.994463 

224.60 

1.82 

Keck 

5437.132373 

241.94 

1.83 

Keck 

5439.129924 

226.52 

1.75 

Keck 

5440.044256 

191.75 

1.83 

Keck 

5442.191536 

203.51 

12.40 

OAO/S 

5444.192803 

198.81 

12.47 

OAO/S 

5445.188846 

187.01 

11.86 

OAO/S 

5455.950750 

124.52 

1.96 

Keck 

5465.039075 

148.60 

1.84 

Keck 

5471.108357 

98.93 

11.72 

OAO/S 

5475.112576 

133.07 

14.14 

OAO/S 

5476.113014 

131.94 

13.30 

OAO/S 

5480.084970 

56.51 

17.92 

OAO/S 

5503.986138 

52.05 

11.18 

OAO/S 

5525.008388 

-1.55 

10.30 

OAO/S 

5548.949071 

2.94 

11.34 

OAO/S 

5556.925781 

-52.35 

15.54 

OAO/S 

5570.892782 

-30.05 

15.01 

OAO/F 

5758.302148 

139.44 

12.77 

OAO/S 

5759.268807 

75.19 

15.06 

OAO/S 

5780.140654 

31.75 

7.99 

Subaru 

5812.211520 

-25.19 

9.26 

OAO/F 

5813.142636 

-15.14 

12.14 

OAO/F 

5814.175882 

18.86 

12.09 

OAO/F 

5816.172657 

-9.86 

11.06 

OAO/F 

5817.137730 

-56.90 

9.44 

OAO/F 

5842.058927 

-49.84 

13.19 

OAO/F 

5843.132379 

-70.46 

10.49 

OAO/F 

5845.097980 

-52.67 

10.81 

OAO/F 

5845.119053 

-48.87 

9.46 

OAO/F 

5846.094836 

-31.20 

10.73 

OAO/F 

5846.116016 

-16.82 

11.20 

OAO/F 

5901.014322 

-59.08 

15.15 

OAO/F 

5901.963414 

-65.42 

11.41 

OAO/F 

5903.960129 

-170.84 

11.45 

OAO/F 

5904.963532 

-157.23 

9.10 

OAO/F 

5905.999945 

-153.12 

10.62 

OAO/F 

5907.960761 

-125.94 

8.91 

OAO/F 
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Table 3 — Continued 


Observation Date 
(JD -2450000) 

RV 

(ms -1 ) 

Uncertainties 

(ms -1 ) 

Observatory 

5909.960182 

-125.27 

11.24 

OAO/F 

5910.968214 

-99.10 

17.90 

OAO/F 

5911.934768 

-76.82 

10.49 

OAO/F 

5912.892384 

-85.77 

7.29 

OAO/F 

5912.913863 

-91.26 

8.30 

OAO/F 

5913.903891 

-120.54 

21.31 

OAO/F 

5914.901327 

-168.36 

7.27 

OAO/F 

5914.921310 

-170.06 

13.28 

OAO/F 

5930.879865 

-159.39 

11.11 

OAO/F 

5931.877692 

-182.01 

8.77 

OAO/F 

5932.889575 

-134.12 

8.41 

OAO/F 

5933.885456 

-99.03 

8.10 

OAO/F 

5934.878720 

-102.03 

8.05 

OAO/F 

5935.888159 

-125.08 

12.19 

OAO/F 

6134.296666 

-117.54 

13.65 

OAO/F 

6136.265378 

-90.49 

14.10 

OAO/F 

6136.286852 

-88.34 

12.79 

OAO/F 

6148.092136 

-78.26 

7.32 

Subaru/IS#l 

6149.066180 

-102.46 

6.81 

Subaru/IS#l 

6164.265986 

-101.72 

10.27 

OAO/F 

6165.246902 

-69.81 

10.73 

OAO/F 

6241.051622 

314.86 

15.36 

OAO/F 

6274.952788 

113.00 

7.67 

OAO/F 

6281.899790 

95.35 

8.57 

OAO/F 

6525.307091 

157.21 

10.77 

OAO/F 
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Table 4 

RVs for HD 67087. 


Observation Date 
(JD -2450000) 

RV 

(ms -1 ) 

Uncertainties 

(ms -1 ) 

Observatory 

3361.119675 

85.64 

11.26 

Subaru 

4079.034701 

70.38 

11.03 

Subaru 

4080.041048 

61.08 

8.87 

Subaru 

4081.107316 

63.50 

10.11 

Subaru 

4523.858969 

-27.12 

10.13 

Subaru 

4587.797606 

-81.32 

9.86 

Subaru 

4744.134303 

-31.62 

10.15 

Subaru 

4787.077253 

13.78 

9.21 

Subaru 

5566.282107 

46.39 

13.39 

OAO/F 

5591.177854 

26.16 

16.90 

OAO/F 

5592.185199 

43.64 

13.45 

OAO/F 

5594.176621 

46.84 

17.35 

OAO/F 

5595.178250 

37.70 

13.91 

OAO/F 

5596.178941 

17.20 

11.44 

OAO/F 

5597.200501 

42.93 

13.36 

OAO/F 

5619.140168 

-7.08 

15.94 

OAO/F 

5649.114480 

-15.63 

12.04 

OAO/F 

5651.125157 

-2.28 

13.26 

OAO/F 

5655.035694 

-21.83 

18.42 

OAO/F 

5667.974503 

-53.30 

21.34 

OAO/F 

5669.008561 

-33.89 

14.46 

OAO/F 

5671.004736 

-17.33 

18.93 

OAO/F 

5672.014501 

-26.10 

12.13 

OAO/F 

5695.003326 

-31.91 

14.75 

OAO/F 

5695.960308 

-56.16 

13.23 

OAO/F 

5699.969525 

-39.14 

14.30 

OAO/F 

5700.984984 

-48.07 

16.30 

OAO/F 

5815.320048 

71.48 

17.03 

OAO/F 

5816.321979 

84.62 

17.37 

OAO/F 

5842.301613 

80.10 

11.41 

OAO/F 

5843.336659 

73.52 

12.64 

OAO/F 

5845.304024 

93.58 

15.53 

OAO/F 

5846.320505 

94.84 

12.42 

OAO/F 

5854.102996 

92.03 

8.96 

Subaru/IS#l 

5904.270617 

67.56 

14.79 

OAO/F 

5911.129521 

78.21 

14.21 

OAO/F 

5932.040732 

30.38 

13.50 

OAO/F 

5965.031544 

5.15 

8.41 

Subaru 

5966.079041 

-22.73 

17.53 

OAO/F 

5967.932216 

0.41 

9.97 

Subaru 

5968.821621 

-7.63 

9.58 

Subaru 

6014.031330 

-44.49 

11.49 

OAO/F 

6015.034135 

-51.38 

12.28 

OAO/F 

6040.960413 

-77.30 

12.72 

OAO/F 

6057.947730 

-86.41 

16.10 

OAO/F 

6241.358500 

41.08 

15.18 

OAO/F 

6275.300959 

29.48 

12.55 

OAO/F 

6554.286423 

61.06 

13.23 

OAO/F 

6758.939170 

-105.20 

11.07 

OAO/F 

6780.966475 

-85.21 

11.18 

OAO/F 

6784.989355 

-115.59 

13.17 

OAO/F 
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Table 5 

Orbital Solutions for the newly detected planetary systems 


Parameter 

HD 1605 b 

HD 1605 c 

HD 1666 b 

HD 67087 b HD 67087 c 

p(d) 

577.9+4® 

2Hlt32 

270.0lg-| 

352.211'g 

237411“ 

T C (JD—2450000) 

3443.3+673 

4758.31“®+ 

3526.31s ® 

154.81 2 ®;® 

qon r+370.6 

OZZ.o_4Q7 i 

e 

0.078l°°| 2 5 

0.098l°;°i 

0 rq+0-03 

U.OO—o 02 

0.l7l°;°l 

0 7 fi + 017 

U - ‘ °-0.24 

a; (rad) 

0.45t°;®l 

4.20l°;l° 

c 1 q+0.09 

4.97!°;°® 

4 4«+ 0 - 48 
^°-0.43 

A'i(m s -1 ) 

19.81J; 2 

46.511® 

199.4l|;| 

73-61 it 

qo 0 + 151.4 

yo.o_ 40 s 

a± sini(10 —3 AU) 

i n^ 0 - 06 

J--UD_0.05 

8.98l°;°® 

q oq + 0.13 

O.OO-0.08 

9 qc:+0-18 

Z.OO—016 

13.41“;® 

a p (AU) 

t-48_o.oi 

q ro+0.05 
O.OZ—0 05 

U.y^_o 02 

1.U8-0 04 

3.86l°;« 

Mpsini(Mjup) 

u.yD_o.o4 

3-48lg;l? 

6 4Q+0.3i 

o.^0-0,22 

Q Ofi”*” 0,22 
9.UO-0 20 

4-85ll° 6 ? 

Trend(m s —1 yr —1 ) 

—4 Q2 +0 - 17 
^• yz -0.22 




Nobs 


92 

99 


51 

ffjitt (m s" 1 ) 


6.0 

34.6 


4.0 

rms(m s —1 ) 


6.4 

35.6 


11.8 


0.97 

1.00 

0.90 a 

7HDS&HIDES ( m s 0 


6.7 

-7.4 


17.0 

7hir.es ( m s_1 ) 


0.1 

54.3 




a The best-fit parameters from the MCMC-fit is applied. 



